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The  scattering  of  normal  incident  time-harmonic  TEM  electromagnetic  wave 

by  a  cylindrical  target  with  axis  along  2-direction  is  considered,  e.g., 

.E  “  El  +  El  and  H  ■  H  1  ,  where  i  is  the  unit  vector  in  the  a-direction. 
x- x  y— y  —  r—2  —a  — 

The  whole  space  domain  12  is  divided  into  three 
connected  but  non-overlapping  sub-domains, 
the  interior  region  12^  representing  the 
target  and  possessing  a  non-ortho gonal 
cylindrical  grid  system  centered  in 
Itself,  the  Intermediate  region  12^ 
representing  the  free  space  just  outside 
of  the  target  and  possessing  the  same 
grid  system,  and  the  exterior  region  12^ 
representing  the  far-field  free  space 
but  truncated  at  a  large  distance  away 
from  the  target  and  possessing  the 
standard  orthogonal  cylindrical  grid 
system  (Fig.  1). 

To  facilitate  the  discretization 
of  the  Maxwell's  equations  on  the  non- 
orthogonal  grid  system,  the  following 
Integral  forms  of  the  Maxwell's  equations  are  used 

£l*  dA  ■  iw  ftuHj  *ds, 

$  H^*  •  ft (<T  -  lui£)  *_E^  *  d£, 

'j>  e_2  *  dA  “  i0J  &  uoi2  * 

e 

$  H2-  d l  -  -iweQ  ft  E2  •  ds, 

£  E^*  dA  "  i<A>Uo  ftilj  *  d— ’ 
i  Hj*  dZ  -  ft  (J  -  iue^)  •  ds, 
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with  boundary  conditions  at  3(2 


n  *  I2  "  -  X  -1  ’ 

eo^2*—  ’  Hi*». 
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n  x  H2  -  n  x  H 


and  tha  asymptotic  terminating  condition  at  3ft^» 


n  x  £  -  -  (u  /e  K(n  x  h  ) , 

— J  0  0“  — j 


where  n  is  the  unit  normal  vector  at  the  Interfaces  and 


e  •  0 


U  -  0 


Eq.  (1)  is  discretized  by  using  the  rectangle  rule  on  the  line  integral 
around  the  edges  of  all  incremental  quadrilateral  defined  by  the  grid  system. 
Lee  each  grid  point  of  tha  non-orthcgonal  polar  grid  system  be  denoted  by 

9j)  ■  (i,J)»  where  "i"  and  "j"  denote  the  i-th  closed  cylindrical  grid 
line  and  the  J-th  radial  grid  line  respectively;  let  the  center  of  the 
quadrilateral  defined  by  (i,J),  (i+l,J),  (i,j+l)  and  (i+l,j+l)  be  denoted  by 
(L+*J»J+*5)  »  let  A g  be  the  incremental  distance  between  the  points  a  and  8, 
and  AAy  be  the  area  of  the  incremental,  quadrilateral  centered  at  y. 

Moreover,  let  i  •  1,2,3, ..., I,  and  i  »  l,2,3,...,J. 
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In  the  neighborhood  of  (i,J)  of  the  typical  discretized  (1)  is 


Ei-*,j(ri+l,j'  ri,j)  ”  Ei-*,j+l(ri+l,j+r  'i.J+l5 

+  ’  Ei, j-Hs^i.j-%  "  ^Uz,i+ii,j44j  h**i,iW*±**i.S +V 

CH  -  H  r-»  J4sin%(6  -  0  ) 

^i-lj.j-*  J-'* - l±i - 3 - 

i+l,j+l  ri+l,j"  ri-l,J+r  ri-l,j 

■  <,l.J+%t“l,J-tJjEi,J+JS  "  lt(El-%,J+  E1+*,.J+  Ei-'s,j+l+  El-H5,}+lKrl,]+l'ri,j)c0a>s(e 

*  l‘<El-1i.J+  Ei+%,3+  *l->i,J+l+  Elt^,j+1) 

'  <Ex,l,j+>»(rl,  J+l*ln  9j+l"  rU  j,in  V“in  9j-% 

tEj,i,rt(tu+iM,V  ri.jeo*Ve“Vi)}  • 

in  the  neighborhood  of  (i,j)  of  the  typical  discretized  (1)  is 

Ei+JS,j(ri+l,j“  ri,j)  ~  Ei+*S» j+1 ( Ci+1  ♦  J+l”  Ti-,1+L> 

+  Ei+1,  j+^i+l.J**  ”  Ei,J4l*Ati#j-H|  " 

<Hi-%,  j+*l  "  Hi+*S,jW  4<ri.,1+l*  ri.1)slnia(91-H~  V 

ri+l,j+l+  ri+l,j_  ri-l,J+r  ri-l,j 

-  -iu)eo{Ai1>j+ljEijj+i5  -  %<r1J+l-  rtJ)  cos%(ej+r  6^ 

Ei+*s,J+  Ei-1s,j+l+  Ei-^,j+l)}; 

and  in  the  neighborhood  of  (i,j)  of  fiy  the  typical  discretized  (1)  is 

Ei+4s.j(ri+i”  ri)  “  Ei+*i,j+i(ri+r  riJ  +  Vi.j^Vr  9j)ri+l 

“  Ei,rt(Vi'  9j)ri  “  ia>uoJs(0j+r  9j)(ri+r  ri)(ri+r  ri)Hi+*i,j+»5, 

^i-^.J+^S  “  Hi-^,jW  4ri  8in^^91.K-  9.j) 

ri+l  “  ri-l 

■  -  iuc  r. (8  8  )E,  . ,i  +  source  terms  due  to  the  presence  of  J, 

o  l  j-fi  j  — 

*h're  “i.rt  '  ri,jH9j+i-  V 

“d  M1-U!,1+<s  *  I,<91+1"  9j)t(ri+l,j"  rl,i)rl,1+l  +  <ri+I,j+r  ri,J+l)rl+I,1*' 


In  this  way,  the  most  natural  finite  difference  discratizatlon  of  the  Maxwell's 
equations  on  a  staggered  grid  system  (Fig.  2)  is  obtained. 


If  the  differences  of  the  material  properties  spread  linearly  across  a 
grid  zone  Instead  of  across  the  interface,  then  there  is  no  need  to  impose  the 
boundary  conditions  (2)  at  the  material  interface,  because  the  boundary  condi¬ 
tion  for  the  tangential  component  of  £  is  satisfied  automatically  and  the  other 
three  boundary  conditions  are  also  satisfied  automatically  but  approximately. 

In  this  way,  there  is  no  cumbersome  instruction  and  treatment  at  the  interface 
to  Slow  down  the  calculation  on  the  computer.  The  discretization  of  the 
terminating  condition  (3)  is 

“  (M0/E0)5s  HI-*s,j+V 
J  “  . . J— 1  • 


To  organize  the  above  discretized  (l)-(3)  into  a  linear  algebraic  system, 
we  first  decompose  the  complex  electromagnetic  fields  into  their  real  and 
imaginary  parts,  i.e.,  E  ■  IJ*  +  IE#  and  H  ■  H*  +  iH#  .  Then  the  three 
discretized  complex  scalar  field  equations  become  six  real  scalar  field 
equations.  Next,  let  the  components  of  the  unknown  field  vector  X  be  arranged 
cyclic  in  "J"  for  each  half  integer  Incremental  increasing  in  "1",  i.e., 


X  - 


h%,»s’h5,  iv*  ,h4,  j-v  /Ei,«*,Ei,  1% . E*,j-%,/Ei,y  Ei,  i% . Ei,  j-y 1 

. ^EI,H,EI,  1%’  *  ’  *  *EI,  J-%*  ^EI,Ji,EI,  14’  *  *  *  ,EI,  J-l^  » 


let  the  known  source  vector  be 
B  -  (0,0, 0,0 . ,bl,b2’ . ,b  )T’ 

where  the  non-zero  components  b's  depend  upon  the  spatial  distribution  of 
the  source  J, 

and  the  system  matrix  A  is  a  non-symmetrlc  band-structured  sparse  matrix 
which  is  shown  in  the  following: 

Let  be  defined  by  i  ■  1,2,3,...,!,  defined  by  i  •  T,T+l,...,y, 

and  defined  by  i  ■  y,y+l,  ...,I. 


where  sparse 

“‘ricMi,i  ■  Ai,2  ■  •••  "  4l,T  corr'!,po°‘i  t0  (1>J>  l"  ni-  42.T+1  •  h,x+2 
- -  =2 ,Y  “trMponi  to  (1.J)  In  Bj.  and  ^  ^  -  ...  - 

l4  Aj  j  correspond  to  (i.j)  in 

Because  AX*  li  will  be  solved  many  times  for  different  B/ s,  the  method 
of  LU-decooposition  will  be  used  to  solve  this  linear  algebraic  system  for 
saving  computer  times.  At  this  moment,  serious  programming  effort  has  just 
begun. 


